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Abstract 

We solve the inverse scattering problem for multidimensional vec- 
tor fields and we use this result to construct the formal solution of 
the Cauchy problem for the second heavenly equation of Plebanski, 
a scalar nonlinear partial differential equation in four dimensions rel- 
evant in General Relativity, which arises from the commutation of 
multidimensional Hamiltonian vector fields. 

1 Introduction 

In this paper we solve the inverse scattering problem for multidimensional 
vector fields and we use this result to construct the formal solution of the 
Cauchy problem for the real second heavenly equation 

o,tx-o, zy + e tXX e tyy -e 2 xy = o, e = e(x, y ,z,t)eR, x, y ,z,teR, (1) 

where 8 :X = 89/ dx, 9 >xy = d 2 6/dxdy. 

This scalar nonlinear partial differential equation (PDE) in 4 independent 
variables x, y, z, t, introduced in P by Plebanski, describes the Einstein field 
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equations that govern self-dual gravitational fields. Its 2-dimensional reduc- 
tion 9 jZ = 1 = is the Monge-Ampere equation, relevant in Differential 
Geometry. As we shall see in the following, the heavenly equation plays also 
a distinguished role in the theory of commuting 2-dimensional, Hamiltonian 
dynamical systems. 

The heavenly equation, together with the equations for the self-dual (and 
ant i- self- dual) Yang-Mills (SDYM) fields [2, are perhaps the most distin- 
guished examples of nonlinear PDEs in more than three independent vari- 
ables arising as commutativity conditions of linear operators ||] , and there- 
fore amenable, in principle, to exact treatments based on the spectral theory 
of those operators If the SDYM equations are considered on an ab- 

stract Lie algebra, then the heavenly equation can actually be interpreted 
as a distinguished realization of the SDYM equations, corresponding to the 
Lie algebra of divergence free vector fields independent of the SDYM coordi- 
nates jSJ. Equation ((TJ has been investigated within the twistor approach in 
[7j. A bi-Hamiltonian formulation, a hodograph transformation, and a non- 
linear dressing formalism for equation (|TJ) have been recently constructed, 
respectively, in [S], and [TU] . 

We end this introductory section remarking that the inverse scattering 
problem for vector fields developed in this paper presents some interesting 
new features with respect to the traditional Inverse Scattering Transform 
(1ST) in multidimensions. For instance: 

• The space of eigenfunctions of the multidimensional vector fields is a 
ring. In the heavenly reduction, the space of eigenfunctions is not only 
a ring, but also a Lie algebra, with Lie bracket given by the natural 
Poisson bracket. 

• Constructing solutions of the integrable multidimensional PDEs arising 
from the commutation of vector fields is equivalent to characterizing 
commuting pairs of dynamical systems, Hamiltonian in the heavenly 
reduction. 

• The scattering data are intimately related to the scattering vector A, 
associated with the dynamical system which describes the character- 
istic curves of the vector field. These scattering data, in one to one 
correspondence with the coefficients of the vector field, depend on a 
number of spectral variables equal to the number of space variables. 
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Thus the counting is consistent (a highly non generic feature in multi- 
dimensional 1ST). 

• A central role in this 1ST is played by a certain Riemann-Hilbert prob- 
lem with a shift depending on the scattering vector A. 

• The above ring property allows one to construct two essentially dif- 
ferent 1ST schemes. The first one, in terms of a set of eigenf unctions 
intimately related to the integral curves of the associated vector field, is 
characterized by a nonlinear inverse problem. The second one, involv- 
ing more traditional Jost and analytic eigenf unctions, obtained "ex- 
ponentiating" the above eigenf unctions, is instead characterized by a 
linear inverse problem. 



2 The Cauchy problem 

A multidimensional Lax pair. It is known that the commutation of 
multidimensional vector fields leads to nonlinear first order multidimensional 
PDEs (see, f.L, 0). 

Consider, for instance, the pair of operators 

N 

Li = d t i + \d z i + Ui d x k = d t i + \d z i + Ui ■ Vj, i — 1, 2 (2) 

k=l 

where d x denotes partial differentiation with respect to the generic variable x, 
x = (x 1 , .., x N ), = (d x i, .., 8 x n), Ui = (uj, --,uf), i — 1, 2, A is a complex 
parameter and the vector coefficients Ui depend on the independent variables 
f,z\x k , i = 1,2, k = 1,..,N, but not on A. The existence of a common 
eigenfunction / for the operators L\ and L 2 : 

Uf = L 2 f = 0, (3) 

implies their commutation, VA: 

[L 1 ,L 2 )=0, (4) 

which is equivalent to the following system of 2N first order quasi-linear 
PDEs in (4 + N) dimensions: 

u i,z*= u 2, z i> k = l,..,N, 

U lfi - U 2t^ + E ( U 2 U 1 X 1 - U [ U 2 xi) = °> k = l,..,N. ( ' J) 

1=1 V ' ' ' 
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Parametrizing the first set of equations in terms of the potentials U k 

< = U% i = l,2, k = l,..,N, (6) 

one obtains the following system of N nonlinear PDEs for the N dependent 
variables U k in (4 + N) dimensions: 

N 

Ufa - Ufa + £ {U^Ufa - U^Ufa) = 0, k = 1, .., N. (7) 



i=i 



This system admits a natural reduction; indeed, applying the operator Ylk=i Q 
to equations (J7J) one obtains 



■ k 



N 



d t id z 2 - d t 2d z i + {U k ^d z 2d xk - U k jZ 2d z id xk ) 



k=l 

from which one infers that the condition 

N 



N 



E U U = °> (8) 



1=1 



E^ = ° ( 9 ) 



fc=l 

is an admissible reduction for equation (J7j), implying that the condition of 
zero-divergence : 

V$-Ui=Q, i = l,2 (10) 

is an admissible constraint for the vectors Ui, i = 1, 2. 

From now on, we concentrate our attention on the following important 
example: 

AT = 2, z l = x\ i = l,2. (11) 

In this case, the zero-divergence reduction (jlUj) makes the two vector fields 
Ui-Vg Hamiltonian, allowing for the introduction of two Hamiltonians Hi, i = 
1, 2 such that: 

»; "'7/,.,,, i,3,k = 1,2, (12) 

where eP is the totally anti-symmetric tensor, which, due to (JUt), are parametrized 
by a single potential 9: 

II, 0,- K ' ,A// -- < !>: 0,r- U < jh Hj <' J 0,>. i,j = 1,2). 

(13) 
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Then the compatible linear problems © can be written down as Hamilton 
equations with respect to the times t 1 ,^ 2 : 

f >t i = {H x + \x 2 ,f} 3 , , . 

f ) * = {H 2 -Xx 1 ,f} 3 , U4j 

where {-, - }g is the Poisson bracket with respect to the variables x , x 2 : 

{f,g}x = f, x ig,x2 - f,xig,xi, (15) 

and the nonlinear system (JJJ) reduces to the heavenly equation in Hamiltonian 
form 

6 )t 2 x i — Q )t i x 2 + {9 tX i, 8 iX s}s = constant, (16) 

equivalent to (JTJ) after choosing the constant to be zero, with the following 
identification of the variables 

t 1 = z, t 2 = t, x 1 = x, x 2 = y, (17) 

which we are going to make in the rest of the paper. 

Commuting dynamical systems. It is well-known (see, f.i., |llj ) that 
linear first order PDEs like Q are intimately related to systems of ordinary 
differential equations describing their characteristic curves. The dynamical 
systems associated with the vector fields Li 2 : 

Li = d z + Xd x + ui ■ Vaf, / lg s 
L > 0, ■ Xd y + u 2 ■ V, 



where x = (x, y) and = (d x , d y ), are: 



U : f = «i(f,z) + (A,0), 



U : f = « 2 (x,t) + (0,A). 



(19) 



We remark that the two flows generated by the times z and t commute, 
VA, iff the fields u\ satisfy the integrable quasi-linear equations (|5|L(jrT|). 
Therefore: 

Constructing solutions of the integrable PDEs /f3|)./ lii)) is equivalent to solv- 
ing the classical problem of constructing pairs of commuting 2- dimensional 
dynamical systems. 
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We also remark that, in the heavenly zero- divergence reduction (jlUj) . the 
commuting flows (|19|) are Hamiltonian: 

§ = {x,H 1 + \y} s , 

(20) 

f = {x,H 2 -\x} s . 

Consider now a solution (r 1 , r 2 ) of equation (fTTIk ) and assume that u\ —* 
as \z\ — > oo; then the phase space point f — (r 1 ,?- 2 ) travels asymptotically 
with constant speed (A,0): 

r^s± + (\,0)z, z->±oo, (21) 

where s± are constant vectors. 

If si is given, then s+ can be viewed as a function of si and A: s + (sl, A), 
and the difference A between these asymptotic positions 

A (si, A) = s+ - si (22) 

describes the z-scattering of the phase space point of the dynamical sys- 
tem (|19b). As we shall see in the following, the scattering vector A is also 
intimately related to the z-scattering datum S associated with the partial 
differential operator L\ (see (I51jl ). which plays a central role in the 1ST for 
the vector field L\, and, consequently, in the solution of the Cauchy problem 
for the nonlinear PDEs (|5 |) .([li p . and for the heavenly reduction (jTfl. 

Now we consider such a Cauchy problem, within the class of rapidly 
decreasing real potentials u\: 

u\ -> 0, (x 2 + y 2 + z 2 ) -»• oo, , > 

u{eR, (x,y,z)eR 3 , t>0, 1 ' 

interpreting t as time and the other three variables x, y, z as space variables. 

To solve such a Cauchy problem by the 1ST method (see, f.i., |1],[S]), we 
construct the 1ST for the operator L\ in (JTSk ). within the class of rapidly 
decreasing real potentials, interpreting the operator L 2 in ()18b) as the time 
operator. 

Eigenfunctions and their properties. Since are linear, first order, 
partial differential operators with scalar coefficients, the product of two so- 
lutions /i )2 of the equation Z^/ = is also a solution: L i (f 1 f 2 ) = (in 
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general, an arbitrary function G(fi, fa) of these two solutions is also a so- 
lution: LiG(fi, f 2 ) = 0). Therefore the space of solutions of the equation 
Lif = forms a ring. As we shall see in the following, this property intro- 
duces important novelties in the Inverse Scattering formalism. 

The localization of the vector potential U\ implies that, if / is a 
solution of Lxf — 0, then 

f(x,z, A) -»• f±(X, A), z -> ±oo, . . 

X := x - (A,0)z; 1 J 

i.e., asymptotically, / is an arbitrary function of (x — Xz), y and A. 

A central role in the theory is played by the real vector eigenfunctions 
tp±(x, z, A), the solutions of Li<p± = defined by the asymptotics 

(p±(x, z, A) — > X, z — >• ±oo. (25) 

Their connection to the dynamical system ()19h ) is immediate. Let x = 
r±(u, z, A) be the solutions of the dynamical system (fT§k) satisfying: 

$=«i(r±,^+(A,0), (2g) 
r±(d;, z, A) — > uj + (A, 0)z, z — > ±oo, 

where w 6 R 2 is a given constant. Then the vector eigenfunctions arise 
from solving the equation x = r±(u, z, A) with respect to uj: 

x — r±{ti, z, A) -v=> = <p±(x, z, A). (27) 

Hereafter, when systems of equations are inverted, we will always consider 
a neighborhood of a point in which the corresponding Jacobian does not 
vanish. 

Due to the ring property of the space of eigenfunctions, an arbitrary 
function of 0± is also an eigenfunction. The Jost eigenfunctions (p^ are 
defined by: 

± (f, z- a, A) := e**fe(*.*.A). (28) 
They satisfy the asymptotics: 

± (f,^;a,A) -»• e lS -*, z -»• ±oo, (29) 
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where a = (a>i, a 2 ) are arbitrary real parameters and a- X is the usual scalar 
product in M?. They are equivalently characterized by the integral equations 



p iaX 
e 1 



^(x, z; a, A)+ 

L 3 dx'dz'G ± (x — x',z — z'\ A) (ui(x', z') ■ V^) ± (x', 2;'; a, A) 

(30) 

in terms of the Green's functions 

G ± {x, z- A) = T 0(^z)5(x - Xz)5{y). (31) 

A crucial role in inverse scattering is also played by analytic eigenfunc- 
tions. In the 1ST for the vector field L l5 the analytic eigenfunctions ip±(x, z\ a, A) 
are the solutions of L\^)± = satisfying the integral equations 

ip±(x, z; a, A)+ 

J R3 dx'dz'G±(x -x',z- z'\ \)e id ^-^ (u^x 1 , z') ■ SJ S >) ip±(x', z'; a, A) = e ia ^ , 

(32) 

where G± are the analytic Green's functions 

G± < g ^ ;A) = ± M»-(A ±fe M - (33) 

The analyticity properties of G±(x,z; A) in the complex A-plane imply that 
ip+(x, z; a, \)e~ ia ' x and tf}_( A)e * are analytic, respectively, in the 

upper and lower half A - plane, with the following asymptotics, for large A: 

ij ± (x,z;a,X)e- i ^ = l + ^^^ + 0(\- 2 ), |A| » 1, (34) 

A 

where: 

Q±(x,z) = ±P [ X Z u x {x',y,z') - ^ I / - / I dx'ui(x',y, z), 




; -Mi X , V, Z 

V»o x / 

(35) 

entailing that 

ux(x, z) = iQ± x (x, z). (36) 

It is important to remark that the analytic Green's functions ()33p exhibit 
the following asymptotics for z — > ±00: 



G+(x, z: A) — > ± „ ., S( ~ y ) — n , z — > +00 



G±(£,z;A) ^± 27Ti[x lt±ie] i 
8 



(37) 



entailing that the z = +00 asymptotics of ip + and are analytic respectively 
in the lower and upper halves of the complex plane [x—Xz), while the z = —00 
asymptotics of ip + and ip- are analytic respectively in the upper and lower 
halves of the complex plane (x — Xz) (similar features were obtained in |12j). 

Since the solution space of the equation L\f = is a ring, then the 
product of two Jost (analytic) solutions corresponding to different vector 
parameters a, ft G IR 2 is still a Jost (analytic) solution satisfying: 

± (f, z; a, A)0 ± (f, z; ft, A) = ± (f , z; a + /3, A), . gg . 
ip±{x, z; a, X)ip±(x, z; ft, A) = if)±(x, z; a + ft, A). 

Choosing a = (1,0) and ft = (0,1), equations (JHHj) imply that the Jost 
(analytic) eigenfunctions depend on a only at the exponents: 



0±(f, z; a, A) = (<j>f{x, z; A))" 1 (0±(f, z; A))" 2 = e ^±(^) j ^ 
il>±(x, z; a, A) = (ip ±1 (x, z- X)) ai {ip± 2 (x, z; A))" 2 =: e & ^ 3 ^ ' 



where ip±i(x, z; A) and i()± 2 (x, z; A) are the analytic eigenfunctions satisfying 
the integral equations for a = (1, 0) and a = (0, 1) respectively. 

Due to the ring property, the vector functions jr±, appearing in ()39b) and 
defined by 

tt±(x,z,X) := ( \ogip ±1 (x, z; A), log i>± 2 (x, z; A)), (40) 

are also eigenfunctions of Li : Li7r± = 0. In addition, since if>± are analytic in 
A, V a G M 2 , it follows that t/>i,2 in (139b ) cannot have zeroes and poles in their 
analyticity domains. Therefore tt + and 7r_ are also analytic, respectively, in 
the upper and lower halves of the complex A-plane, with asymptotics: 

tt ± (x, z, A) = iX + Q± (*' Z) + 0(A- 2 ). |A|»1, (41) 

A 

Analogously, the z = +00 asymptotics of 7r + and 7?_ are analytic respectively 
in the lower and upper halves of the complex plane {x—Xz), while the z = —00 
asymptotics of 7? + and 7?_ are analytic respectively in the upper and lower 
halves of the complex plane (x — Xz). 

Since the vector fields are real (ui = u*), we have the following obvious 
restrictions for the solutions ^ and ip± (for real A): 

</> ± *(a;, z; a, A) = ± (x, z; — a, A), (42) 
ip±*(x, z; a, A) = ip T (x, z; — a, A). (43) 
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As we shall see in the following sections, the eigenf unctions (p± and tt± are 
the basic ingredients of an 1ST involving a nonlinear inverse problem, while 
the eigenf unctions ± and ip± are the basic ingredients of an 1ST involving 
a linear inverse problem. 

Scattering data. We begin this section observing that the vectors f± de- 
fined in (|26p. can be identified with the vector r in (|21|) setting s± = u; then 
the asymptotics of f± at the opposite ends of the z axis can be described in 
terms of the scattering vector A in (|22|) as follows: 

r-(u,z, A) -> s+(uJ, A) + (A, 0)z, 2^+00, 
f+(u;, z, A) — > 0(cJ, A) + (A, 0)z, z — ► —00, 

where 

s + (u, A) = u5 + A(aJ, A), (45) 
and a> = fi(s+, A) is the inverse of the transformation from cJ to A): 

s + = uj + A(u, A) <^> u = tt(s + ,\). (46) 
Due to the intimate connection (|27j) between r ± and </?±, also the asymptotics 

— * 

of (p± at the opposite ends of the z axis are described in terms of A: 

<P-(x, z, A) -> fi(A, A), ^ -> +00, 
<p+(x, z, A) — > s + (A, A), z — » — 00, 

and, consequently, the expression of y? + in terms of y?_ is given by: 

(p + (x, z, A) = s + (ip-(x, z, A), A) = y?_(f, z, A) + A(<£_(x, z, A), A). (48) 

To derive equations (|47|). it is sufficient to take the 2; = ±00 limit of (J27J), 
using (|44|) : equation (|4*Hj) follows from the ring property and from ([47)1 . 

The Jost solutions form a natural basis for expanding any eigenfunction 
of L\. The expansion of (f) + in terms of <fi~: 

<f) + {x,z]a,\) = / d0S (a, P, \)<f>~(x, z; P, A) (49) 

defines the ^-scattering datum S 1 , which admits the standard integral repre- 
sentation: 

5(5, (3, A) = - P) + / dxdze-^* {u^x, z) ■ V*) + (f , 2; a, A) (50) 
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following from (fSUJ). where 5(a — (3) = 5(ot\ — /3i)5(a 2 — 

The ^-scattering datum S exhibits the following Fourier representation 
in terms of the scattering vector A: 

S(Sj,X)= [ J^ e ^-^^), (51) 

This basic formula is consequence of (J2EJ), (}4"5j) and pOjl. Indeed, substituting 
d2Hl± hito one obtains: 

e***+ = g(3 )V >_, A), (52) 
where 5 is the Fourier transform of 5, with respect to (3: 

S(a,u,\)= [ d$S(aJ,\)(*fa. (53) 

Then, using (|48|h it follows that 5(a, a;, A) = e 4a - s +(^> A ). Therefore the special 
dependence (|2*%|) of on a only at the exponent implies that also the Fourier 
transform S of S depends on a only at the exponent. 

The direct problem is the mapping, via (|5U|). from the two real potentials 
u\ 2 , functions of the three real variables (x,z), to the scattering datum S, 
or, via (|51jl. to the two real components A 1 ' 2 of A, functions of the three 
real variables (u, A). Then, thanks to the ring property of the space of 
eigenfunctions, the counting is consistent. 

The expansions of the analytic eigenfunctions ip± in terms of the Jost 
eigenfunctions: 

ip + (x, z; a, A) = J R2 df3K±(a, (3, A)0 ± (f, z; (3, A), ^ 
■0_ (if, 2; a, A) = L 2 dpK^(a, (3, A)0 ± (x, 2; /3, A) 

introduce the kernels K± (a, (3, A) ; they are connected to 5 via the integral 
equations 

K±(aJ,\)= [ d1K±(an,X)StfJ,X), (55) 

which follow directly from (|54j) . replacing <p + by its expansion f|49|) in terms 
of</r. 

The reality of the vector fields and its consequence (|42|) imply the follow- 
ing reality constraints for the scattering data S and K (for real A): 

S*(a, (3, A) = S(—a, -(3, A), (56) 
K±*(a,P,\)=K±(-a,-j3,\). (57) 



11 



Due to the special dependence f)39p of if)± and <p± on a, also the Fourier 
transforms K±(a, to, A) of K±(a, (3, A) with respect to /3, depend on a only 
at the exponents: 

Ki(a,u,\)=(Ki\ujA)) (Ki\u,\)) , (58) 
yielding the following Fourier representations of the kernels K 

K±(aJ,\) = J R2 *fr<&<*-ft 



(59) 



where 

Xl = {xl\xf), xV = ^K^ -iu ]} j = 1,2, (60) 

and yielding also the representations of the analytic eigenfunctions tv± in 
terms of <p± and x± : 



7f+(f, z, A) = i0±(x, z, A) + x+(^±(x, z, A), A) 
7f_(£, A) = i<p±(x, z, A) + x_(<p±(£, A), A) 



The proof is the same as before; replacing the exponential representations 
()28|) and P9*b) of ^ and -?/>±, f.i., in (153k ). one obtains the relation 

e s -*+ = K±(a,<p±,\), (62) 

proving that the dependence of K+ on a is only at the exponent. In addition, 
using the expression of K± in terms of x+i one obtains (|UTk ). 

We have clarified the impact of the ring property on the kernels K± . The 
analyticity properties of the z = ±oo asymptotics of ifi±, discussed in the 
previous section, imply instead the following triangular structures: 

K±(a, 0, A) = 6(a - 0) + 9(±(a 1 - (3 x ))K±{aJ, A), 
K±(a,0,X) = 6(a-/3) + 0(=f(c*i - Pi))K±(aJ, X) 1 ' 

(the same features were observed in ^2])- Therefore the linear integral equa- 
tion ({55} describes a factorization problem allowing to construct, in principle, 
the kernels K± in terms of the given scattering datum S. 

Comparing (|53j) and ()58|) . one infers that KZ, KZ are analytic in the 
upper half u 1 - plane and Kf, KZ are analytic in the lower half uj 1 - plane. 
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Since these analyticity properties are valid V a G R 2 , it follows that K^ 1 , K^ 2 
in (|SB} cannot have poles and zeroes in their analyticity domains. Therefore 
also x± ' 2 (a3, A) an d x± 1,2 (^)^) are analytic, rispectively, in the upper and 
lower half uu% - plane, with asymptotics: 

X | 1,2 (a5,A)^0, \u\ 2 ^oo, A^O. (64) 

At last, plugging the Fourier representations (|BT|) and (|59p into the integral 
equation f|55j) . one obtains the following basic equations: 

Xt(c0, A) - xT(o; + A(u), A), A) = iA^(uj, A), 

£Z(aJ, A) - x-{£ + A), A) = zA 1 ' 2 ^, A), e M 2 , AeR 2 , Xe R, 

(65) 

wich must be viewed as "inhomogeneous Riemann - Hilbert problems on the 
line Im uj 1 = 0, with the shift A 1 (o;, A)" ^3JE] (the real variables uj 2 and A 
appear just as parameters). Given A(uj, A), namely, given the shift and the 
inhomogenuity in (|B3)l. one can construct, in principle, the unique solutions 
X± of these Riemann - Hilbert problems, expressed in terms of the following 
linear integral equations of Fredholm type : 

Xtiu 1 , ^ 2 , A) + JL / R drM{uj\ r, A)x+(r, uj 2 , A) + L(uj\ uj 2 , A) = 0, 
T-{uj\ uj 2 , A) - i / R rfrM^ 1 , r, A)xl(r, ^ 2 , A) + L^ 1 , ^ 2 , A) = 0, 1 ) 

where 

ar/, ,1 , , ,2 n _ dsl(T,u 2 X)/dr 1 

L(uj\ uj 2 , A) = uj 2 , A) + / R dr ^^^l^ Air, uj 2 , A) 

s^w 1 , w 2 , A) = uj 1 + A\uj 1 , uj 2 , A). 

(67) 

The inverse problem. Once x± are known, one reconstructs the eigen- 
function y?_ solving the following nonlinear integral equation 

A) + PiX+{0_(£, z, A), A) + P+xZ{0-(x, z, A), A) = iX, (68) 

where are respectively the (+) and (— ) analyticity projectors with respect 
to A: 

- + 1 f dX' 
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Equation (|68|) is obtained applying P x and P x respectively to the analytic 
(and decreasing at A = oo) expressions (n + —iX) and (ff_ —iX), subtracting 
the resulting equations, and using (JBTJ). 

From the solution of (|68|h one reconstruct the potentials U\ using the 
formulas 

U!(x , z) = id x lim ( A[z(<^_(f, z, A) - X) + x+ (<£_(£, z, A), A)] ), (70) 

consequence of (JSfijl and (jUJ. 

Therefore we have constructed an 1ST which involves a nonlinear step: 
the nonlinear equation (jUHj) of the inverse problem. The above results can be 
summarized as follows. 

Nonlinear 1ST scheme. The direct scattering is the mapping from the po- 
tential ui(x, z) to the scattering vector A(u, A). Knowing A(u, A), one con- 
structs the vectors x+(o>, A) and xZ(u, A) solving the linear RH problems with 
shift In the inverse problem, from the knowledge of x^(o3, A), xl(o>, 
one solves the nonlinear integral equation for <p-(x, z, A) and, finally, 
one reconstructs the potentials U\{x,z) from (|7Uj). 

It is quite remarkable that the inverse problem becomes linear, if ex- 
pressed in terms of the Jost and analytic eigenf unctions ± and ip±. Indeed, 
from the solutions x± °f the R-H problems (|55Jl one can construct, via (|55)l. 
the factorization kernels Xj:. Known KZ,KZ, one reconstructs the Jost 
eigenfunction 0~ by solving the following linear integral equation: 

0-(cf, A)e-^ + P+ J R2 dp6{-{ai - p 1 ))KZ(Sj 1 \)4r0, X)e' &ji + 
PZ J R2 dp9(-( ai -/3 1 ))^(a,/3,A)r(/3,A)e-^ = 1. 

(71) 

This equation, in which we have omitted, for simplicity, the parametric de- 
pendence on (x, z), is consequence of the analyticity properties of ip±, and is 
obtained multiplying the equations 

ip + (x, z; a, A) = f R2 d(3KZ(a, (3, A)(j)~(x, z; (3, A), ^ 
i/}-(x, z; a, A) = / R2 d/3KZ(a, (3, \)<j)~{x, z; (3, A) 

in by e~ 4a ' x , subtracting 1, applying respectively P~ and P^, and adding 
the resulting equations. This inversion procedure has been already presented 
in 
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Once ± are known and, via (J72j) . ip are also known, the potential U\ is 
reconstructed, via (J53J) and ([3T)j). through the formulae: 



2;) 



ic^a; lim A (V'i 6 



A— >oo 



i0 a lim A ( i)±e~ id - x - 1 

A^oo 



- 1 



cf=(l,0) 
a=(0,l) 



(73) 



This second version of the 1ST consists of linear steps only, and can be 
summarized as follows. 

Linear 1ST scheme. The direct scattering is the mapping from the potential 
Ui(x, z) to the scattering vector A(u3, A). Knowing A (63, A), one constructs 
the vectors x+{uj, A) and x~_(c<3, A) solving the linear RH problems with shift 
(Ifi5j) . Then, using the Fourier representations ()59j) . one obtains the kernels 
iTZ and K+. In the inverse problem, the Jost eigenfunction is then recon- 
structed through the linear integral equation (J7T|) , the analytic eigenfunctions 
via (J72J), and the potentials via (J7HJ). 

The small field limit. In the small field limit \u{ ,2 \ = 0(e) << 1, the 
scattering vector A is expressed in terms of u\ in the following way: 

A(uj,A)= f dzu^cu 1 + \z,co 2 ,z) + 0(e 2 ), (74) 
Jr 

while the Riemann - Hilbert problems with shift (J65|) reduce to the scalar 
Riemann - Hilbert problems: 

Xt(u, A) - x^(u, A) = xl(", A) - £t(tf, A) = *A(cJ, A) + 0(e 2 ), (75) 

entailing 



A) = i J R2 pgft^ + A*, a*, z) + 0(e 2 ), 
X^, A) = -£ / Ra jzggjaytfitt + Az, c 2 , z) + 0(e 2 



(76) 



At last, the small field limit of the inverse problem equations (17 lj) and (|f)8j) 

reads 



•(</?_(£, *, A) - X) + P+xl(X, A) + P^xX(X, A) = 0(e 2 ). (77) 
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The t-evolution. The time evolution is introduced in the usual way. We 
observe that the Jost eigenfunction + solves the equation L 2 <f) + = i\a 2 (p + 
which, evaluated at z — — oo using also (JH5J), leads to 

S(a, j3, A, t) = S{a, (3, A, O^ 03- ** (78) 

From (J78J) and ()55|) one infers the same elementary t-dependence for K±: 

K±{aJ,\,t) = K±{aJ,\,0)e iX{a2 - Mt . (79) 

At last, using (fTSj) and the Fourier representation (joTj) of S, one obtains the 
t-evolution of the scattering vector A: 

A(u,\,t) = A(w- (0,A)t,A). (80) 



The heavenly reduction. The heavenly equation corresponds to the Hamil- 
tonian reduction (JTUJ). In this case, if /i 2 are solutions of I/j/i,2 = 0, then the 
Poisson bracket (fTB} of /i j2 is also a solution: Li({fi, f 2 }$) = 0. Therefore 
the solution space of the equation Lif = is not only a ring, but also a 
Lie algebra, with Lie bracket given by the Poisson bracket (jl5|) . This result 
follows immediately from equation (pHk ) and from the Jacobi identity. 
Under the hamiltonian reduction (fTU|) . equation gets the form 

(IT 

— = {r,H 1 + \r 2 },. (81) 

Thus the ^-evolution of r*is canonical, i.e. it preserves the Poisson brackets of 
ri(z) and r 2 (z). So, comparing eqs. lffilb ) and (JUJ), one obtains the following 
constraint for the scattering vector A: 



d(sl,sl) 



1. (82) 



3 Conclusions and open problems 

In this paper we have constructed the 1ST for multidimensional vector fields, 
and we have used it to solve the Cauchy problem for the heavenly equation. 
Although these results have been derived in the particular case (JHJ), directly 
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connected to the heavenly equation, they can be extended in a straightfor- 
ward way to the general vector fields (J2J) and to the (4 + N) dimensional 
PDEs ©. 

Interesting open problems under present investigation are: i) the tran- 
sition from the above formal results to rigorous ones, ii) The study of the 
spectral mechanisms (if any) which could cause a breaking of the initial pro- 
file for the PDEs under investigation (a typical feature of quasi-linear PDEs) . 
iii) The costruction of explicit solutions within this spectral formalism. 
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